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By phenomenologically describing the high-redshift star formation history, i.e., p„(z) oc [(1 + z)/4.5] -a , and semi-analytically 
calculating the fractions of high-redshift Pop I/II and Pop III stars, we investigate the contributions from both high-redshfit Pop I/II 
and Pop III stars to the observed near-infrared (3 yum < A < 5 jim) excess in the cosmic infrared background emission. In order 
to account for the observational level of the near-infrared excess, the power-law index a of the assumed star formation history is 
constrained to within the range of < a < 1 . Such a constraint is obtained under the condition that the virial temperature of dark 
matter halos belongs to the range of 500 K < r vil < 10 4 K. 
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1 Introduction 

In recent years, the cosmic star formation history (CSFH) is 
of more and more interest in both astrophysical and cosmo- 
logical studies, especially the early universe where direct de- 
tections are difficult and scarce. As representative objects 
in the early universe, the first generation stars, commonly 
called Population III (hereafter Pop III) stars, would play a 
key role in the early evolution of the cosmic structure. Pop 
III stars are thought to originate from primordial metal-free 
gas at redshifts exceeding z ~ 10 and precede the normal 
metal-enriched stellar populations (see References [QJ[3] for 
reviews of Pop III stars). Direct observations of Pop III 
stars with current telescopes are impossible due to their high 
redshifts, although these stars are very massive and lumi- 
nous. Alternatively, Pop III stars are expected to have left 
a great amount of diffuse radiation in the universe, which is 
shifted today into the infrared wavelength. For instance, for 
a photon at Lyman-a limit with A s = 912 A emitted from 
z = 20, its wavelength to be observed will be redshifted to 
A r = A s (l+z) = 1.9 |im. Therefore, a significant near-infrared 
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component from Pop III stars is expected to be present in the 
cosmic infrared background (CIB). 

In the aspect of observation, a substantial near-infrared ex- 
cess over the net fluxes produced by galaxies has been in- 
deed inferred from the measurements of CIB anisotropics 
with COBE/DIRBE d, Deep Spitzer Infrared Array Cam- 
era (IRAC) EE), Infrared Telescope in Space (IRTS) Q, and 
Hubble Ultra Deep Field (HUDF) ® . The commonly known 
bolometric flux of such near infrared background excess 
(NIRBE) reads F°^ ess = (2.9 ± 1.3) X 10~ 5 erg s^cnT^r" 1 
||9l . However, a more recent analysis by [8. 10| claimed 
that the flux of the NIRBE within the wavelength range of 
3 fim < A < 5 fim could be ten times smaller as 



^excess ~ 1 X 10 



erg s 



(1) 



because much of the previously estimated excess could be 
due to inaccurate zodiacal light modeling IfTTI . 

Therefore, in this paper, we would use the observed 
NIRBE to constrain the formation history of Pop III stars as 
well as the properties of the first collapsing dark matter ha- 
los hosting such stars. These halos may be too faint to be 
observed with the current telescopes. Different from some 
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previous research studies, in our considerations, we also take 
into account the contribution to the NIRBE from normal Pop 
I/II stars at high redshifts, because in the process of deriva- 
tion of the NIRBE flux, only the emission from the galax- 
ies at relatively low redshifts z < 10 was subtracted (8). In 
the next section, we would present our theoretical consid- 
eration of the formation rates of Pop III and Pop I/II stars, 
where several free parameters are introduced. In section III, 
by confronting the theoretical NIRBE with the observational 
NIRBE, we constrain the model parameters. Finally, a sum- 
mary and discussion are given in the last section. 

2 Star formation rates 

In principle, the formation rate of Pop III stars can be the- 
oretically derived in the hierarchical formation model (e.g., 
|[T2l[T3ll ), which however involves many complicated astro- 
physical issues. Therefore, in order to decrease the model 
uncertainties and utilize some observational information, we 
here give a phenomenological description for the CSFH and 
a semi-analytical calculation for the fraction of Pop III stars. 

2.1 Cosmic star formation history 

Following a series of measurements of star formation rate 
density, especially the complication of Hopkins & Beacom 
fl4"l . a consensus on the CSFH now emerges up to redshift 
z ~ 3.5, which includes a steady increase of star formation 
from z — to z = 1, and a following plateau up to z ~ 3.5. 
An empirical fitting can be written as 



P*(z) 



(l+z) 3M . 
(l + z)°, 



for z < 0.97, 
for 0.97 <z < 3.5, 



(2) 



with a local rate density p*(0) = 0.02 M yr~'Mpc~ 3 . In 
contrast, for higher redshifts, the situation is still ambiguous. 
It could continue to plateau, drop off, or even increase, e.g., 
as considered in |[T"5l . Therefore, we here simply introduce a 
new free parameter a to parameterize the high-redshift CSFH 
as follows 

p t (z) oc {I + zT a , forz>3.5. (3) 

The power-law assumption adopted above is motivated by 
the shapes of the history at z < 3.5 (see equation [2]i and 
also implied by some preliminary measurements of the high- 
redshift star formation rates through Lyman break galax- 
ies lfl6l and through gamma-ray bursts as well as their host 
galaxies HMO). 

As a phenomenological consideration, we further assume 
that the star formation rate density given by equation © has 
contained both the contributions from Pop I/II and Pop III 
stars. 



2.2 Fraction of Pop III stars 

In the hierarchical formation model, star formation takes 
place during the collapse and merging of dark matter halos. 



The formation of the halos hosting zero-metallicity Pop III 
stars should satisfy two fundamental conditions: (i) The ha- 
los must be "freshly" formed through collapse of diffuse dark 
matter but not through merging of smaller halos. In these 
new halos, the stars are formed for the first time, (ii) The 
newly formed halos are not located in the wind radius, R w , 
of old galaxies, because the intergalactic medium within R w 
has been metal enriched by the galactic winds (i.e., feedback 
effect) l2TI - l23l . Here the wind radius R w is defined as the 
Sedov radius of the galactic wind, which is mainly deter- 
mined by the total energy of all supernovae in the galaxy. 
For typical energy release per supernova E$n ~ 10 52 erg and 
typical mass of Pop III stars M ~ 100M G , we can have 
R w ~ 380 (f7z/10 10 M G ) 1/5 [(l + z)/10]- 3/5 kpc. 

For the first condition, following Furlanetto & Loeb 11241 . 
the fraction of matter that can collapse into "new" halos can 
be expressed by 



/new(z) 



f 



n(m min ) 



d 6 c (z') 



dz' o-(m min ) 



dz', (4) 



P cr(m min ) 

where p is the mean cosmic density, 6 c (z) is the critical den- 
sity for collapse, cr(m) 2 is the variance of the density field on 
the mass scale m, n{m) is the Press & Schechter mass func- 
tion 11251 . and m m i n is the minimum halo mass. With a given 
virial temperature T wu , the minimum halo mass can be esti- 
mated by equating the gravitational force of the halos to the 
opposing gas pressure ||26l 

^»»io»«.(^)- w (^f(^r. <» 



10 4 K 



10 



where \i is the mean molecular weight, with \x = 0.6 for ion- 
ized gas and p. = 1.22 for neutral primordial gas. Here we 
choose p - 1.22 as the universe is mostly neutral at high red- 
shifts. 

For the second condition, the probability that a new halo 
lies within R w of an existing galaxy can be approximated 
by JH 



1 



exp 



-(1 + B) 



f 



p 



-r){m)n(m)dm 



(6) 



where T](m) is the ratio of the mass metal-enriched by galactic 
winds to the total mass of each galaxy. Following Furlanetto 
& Loeb 11241 , we have 



T](m) ~ 27 K, 



,1. 



-2/5 



10 



3/5 



(7) 



w \10 10 M o 

where the normalization factor < K w < 1/8 accounts 
for many feedback factors (e.g., radiative losses from su- 
pernova shocks). A higher value of K w represents stronger 
feedback effects. For the formation of Pop III stars, such 
feedback effect is a negative effect. The coefficient B ap- 
pearing in Equation (O gives the excess probability that two 
galaxies sit near each other, by considering that the col- 
lapsing halos are actually biased due to the existing halos 
and inclined to be near the existing halos. To be specific, 
B ~ (mm(Rn)b(m min ){b(m)), where £ mm (R w ) (with R w to be 
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Figure 1 The CSFHs of Pop III stars (dashed lines), Pop I/II stars (dot- 
ted lines), and their sum (solid lines), where a — 1.0 and K w = 0.02 are 
taken. The shaded region represents the variation of the virial temperature 
from T vir = 500 K to 10 4 K. 

the average wind size in the comoving units) is the dark mat- 
ter correlation function 11271 , b(m) is the bias function 1281 , 
and the bracket (• • • ) represents an average on the Press & 
Schechter mass function 



2.3 Formation rates of Pop III and I/II stars 

Combining the two aspects addressed above, the formation 
rate density of Pop III stars as a function of z can be given by 



, df new /dz\ 
Pm(z) = p*(z) — — — (1 - p e ), 
d f con/ dz 



(8) 



where / C oii(z) is the fraction of matter that can collapse into 
"all" of halos. For the Press & Schechter mass function, it 
reads |24) 



/coii(z) = erfc 



(9) 



V^crtmmin) 

Meanwhile, the formation rate density of Pop I/II stars can be 
obtained by subtracting the rate density of Pop III stars from 
the total star formation rate density 



pi/n(z) =P.(z)-pm(z). 



(10) 



In the model described above, three free parameters re- 
main: a, r v ; r , and K w . As theoretically considered, the value 
of TVjj- is in principle determined by the cooling style of the 
halos. The radiative cooling through molecular hydrogen 
would lead to T ylx ~ 500 K B30U31L while the atomic cooling 
gives r v ; r ~ 10 4 K (e.g., I321I33II ). Hence the specific value 
of r v ir relies on the fraction of molecular hydrogen. For the 
temperature range of 500 K < r v ; r < 10 4 K, we present a nu- 
merical result of the CSFH in Figure 1, which indicates that: 
(i) the formation of Pop III stars would dominate the cosmic 
star formation for z > (20 ~ 40), and (ii) higher virial temper- 
ature r v ; r could extend the formation of Pop III stars to lower 
redshifts. 



3 NIRBE emission 

Instead of an elaborate fitting to observations, here we would 
use the observed bolometric flux of the NIRBE (Equation 1) 
to constrain the model parameters. In the following calcula- 
tions, the contributions from Pop III stars and high-redshift 
Pop I/II stars are both taken into account. 



3.1 Pop III stars 

As massive, zero-metallicity stars, Pop III stars would radi- 
ate close to the Eddington limit L e dd = (4nGm p c / 'o~t)M 
1.3 x 1O 38 (M/M ) erg s -1 , where M is the stellar mass. 
Then the lifetime of Pop III stars can be estimated by ti 
eMc 2 /L e( jd - 3.1 x 10 6 yr, which is independent of the stellar 
masses, where s = 0.007 is taken as a fiducial value. Here 
Pop III stars are not considered for significant mass loss dur- 
ing their lifetimes because of their zero metallicity. The black 
body effective temperature of Pop III stars can be estimated 
to Tin ~ 10 5 K with the Eddington luminosity and the theo- 
retical Hertzsprung-Russel diagram of Pop III stars H341I35I . 

At redshift z, the density of the energy radiated from Pop 
III stars can be expressed by 



£m(z) ~ epm(z)c t H (z), 



(ID 



which is released within the duration of ti, where tn(z) ~ 
(2/3)H(z) is the Hubble timescale at that redshift. In the 
standard A-cold dark matter cosmology, the Hubble expan- 
sion rate reads H(z) = Hq y(I + z) 3 Q m + D.^. Hereafter 
H = 71 km s^Mpc" 1 , Q A = 0.73, and Q„, = 0.27 are taken. 
Therefore, the bolometric flux of the NIRBE contributed by 
Pop III stars can be calculated by 



^- I JZm(Z) aV c (z), 



\nd\(zY 



(12) 



where J^m(z) represents the fraction of energy that enters 
the near-infrared wavelength range of 3 pm < A < 5 pm, 
di(z) and dV c (z) are luminosity distance and comoving vol- 
ume element, respectively. To be specific, d[(z) = c(l + 
z) J~ H- l (z')dz' and dV c (z) = 4nd c {z) 2 cH- l {z)dzwithd c (z) = 
di(z)/(l + z). By considering that O a Q m (l + z) 3 for high 
redshifts, we can simplify the above equation to 



cl 



in 



67T//2Q, 



^m(z)pm(z) 
(1+z) 5 



dz, 



(13) 



where Equation (JTTJ is used and lm = Le^a/M = 6.3 x 
10 4 erg s _1 g _1 represents the luminosity per unit stellar mass. 
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Figure 2 The dependence of the correction factor Jt(z) on redshift for 
stellar temperature T = 10 5 K (solid line), which peaks at z p = 24 due to 
Lyman-a- absorption. The case without Lyman-a absorption is also presented 
by the dashed line for comparison. 




Figure 3 The distribution of the ratio F' x ' ccss IF 'x C ' css in the K w — 7\,ir pa- 
rameter space. 



For a black body spectrum described by the Planck func- 
tion Bt(A) for temperature T, the correction factor J?T(z) can 
be generally calculated by 



f 

J A., 



/(!+<-) 



B T (A)dA 



J °°B T (A)dA 



(14) 



with A\ = 3 jt/m and At - 5 fim iflOll . However, as revealed 
by WMAP, the reionization of the universe completes at most 
as early as redshift ~ 10.5, before which the universe is neu- 
tral 11361 . Consequently, for stars at high redshifts, photons 
emitted from them with wavelengths shorter than the Ly-a 
wavelength (i.e., A e < 1216 A) must be all absorbed by the 
diffuse neutral hydrogens (as the Gunn-Peterson trough in the 
spectra of high-redshift quasars). Therefore, the lower limit 
of the integral in the numerator of Equation (fl4l needs to 
be changed to max[/li/(l + z), 1216 A]. In the extreme, if 
<V(1 + z)< 1216 A, then we have JT(z) = 0. 

For Pop III stars with T = 10 s K, a numerical result 
of J£m(z) as a function of redshift is presented in Fig- 
ure [2] As is shown, it firstly increases as (1 + z) 3 up to 
z p = 3 /mi/ 12 16 A - 1 = 24, and then plunges to zero at 
Zmax = 5 /mi/1216 A - 1 = 40. The peak value of J£m(z) at 
z p can be analytically calculated by 



Jftll(Zp) = 



2nhc 2 



(l+z P fAA 



o-T A A 5 c exp[(l +z P )hc/k B TA c ] - 1 



. (15) 



where A c — 4 fim and AA = 2 //m. Moreover, due to 
(1 + z p ) <K 4-ksTAc/hc, the peak value can be approximately 
expressed by 



7ick B AA 1 , 
which is proportional to T~ 3 for T > 2 x 10 4 K. 



(16) 



3.2 Pop I/II stars 

Following Equation ( fT3l , the contribution to the NIRBE 
emission from high-redshift Pop I/II stars can be calculated 
by 



dF, 



i/n 

excess 



cl 



i/n 



dM 



6nH 2 Q Q, 



■I 



(1 + z) 5 dM 



dz, 



(17) 



where M = M/M Q is the dimensionless stellar mass. The 
lower limit of the above integral is set at z ~ 10 since the 
contributions from Pop I/II stars at relatively low redshift 
z < 10 have been subtracted from the excess emission JS). 
The differential form is taken in Equation (fTT) due to the 
wide distribution of the mass of Pop I/II stars. So, differ- 
ent from Pop III stars, both the quantities h/u and J£i/n(z) 
here should be treated as functions of stellar mass. Specifi- 
cally, the empirical mass-luminosity relation for Pop I/II stars 
tells us h/u x M 2 - 5 (L Q /M Q ). Meanwhile, we can estimate 
the stellar temperature by T ~ 5000M 6 K according to 
the Hertzsprung-Russel diagram of Pop I/II stars. For the 
Salpeter initial mass function, <p{M) oc M -2 - 35 , the differen- 
tial star formation rate can be expressed by 



dpywiz) 
dM 



= Ap m (z)M(f>(M), 



where A 



™™ M(p{M)dM 
= 100. Then, we have 



rill 
J.A7„ 

. Tl 

pirn 



0.17 with M n 



(18) 



0.1 and 



Ac 



6nH 2 fl m 



M 



(1+z) 5 



dzdM. (19) 



In Figure 3, we plot the distribution of the ratio of F^ 
to /^excess i n th e K w - T yK parameter space 1 1 which indicates 



m 

excess 



The dependence of this ratio on the parameter a is insensitive, because a influences the formation rates of Pop III and I/II stars in a same way. 
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Figure 4 Dependence of Fexcess on the free model parameters for Pop III stars (solid lines) and high-redshift Pop I/II stars (dashed lines), where the fidu- 
cial values of the parameters are taken to 7" v i,- = 5000 K, K w = 0.01, and a = 0. The horizontal dotted line represents the observed flux of the NIRBE as 
1 X 10~ 6 erg s -1 cirr 2 sr~'. 



that the contribution to the NIRBE from high-redshift Pop I/II 
stars could be much more significant than previously consid- 
ered, especially for low T wu and high K w . 

From Equation (fl9l i, we can see that the value of ^excess 
could be mainly determined by the highest-mass stars having 
temperatures T > 2 x 10 4 K. Therefore, following Equation 
( fToT ), we can get 

■*vn(z) „ (Tw) 3 „ (m 
JtHz) ~\Tm) ~ 8000' 

where the peak redshift z p is replaced by an arbitrary red- 
shift z because the profiles of the functions Jf(z) are basically 
identical so long as T > 2 x 10 4 K. For a more direct com- 
parison with Pop III stars (Equationll3l. we rewrite Equation 

67t^Q m J (1 +Z )5 

where the equivalent luminosity per unit stellar mass reads 

t l \ r m 2 - 95 

l'im=A\ — \ t/M = 850ergs-V, (22) 

\M@)J 8000 s s , K J 

which gives a ratio l'yn/lm = 0.014. Therefore, if pjyn is hun- 
dreds times higher than pm, the contribution to the NIRBE 
from high-redshift Pop I/II stars can exceeds that from Pop 
III stars. As exhibited in Figure 3, such a situation can arise 
with low r vu - and high K w . 

3.3 Constraining the CSFH 

The dependence of the model-predicted NIRBE fluxes on the 
parameters T vn , K w , and a are shown in Figure |4] for both 
Pop III and I/II stars. As can be seen, the uncertainty of the 
NIRBE fluxes mainly comes from the uncertainties of a and 
r v i r , whereas the variation of K w only slightly influences the 
value of F excess . 

Combining Equations (fT3l ) and ( fl9l ), the model-predicted 
flux of the NIRBE emission contributed by both Pop III and 
high-redshift Pop I/II stars can be calculated by 




-0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 5 The parameter region where the model-predicted NIRBE flux 
can account for the observational one (shaded region). The two boundaries 
of the region corresponds to the virial temperatures of 500 K and 10 4 K, 
respectively. 



rmodel plU pL/11 m) 

excess excess excess" V^-V 

As shown in Figure|5] by equating /•'excess to ^excess' we ^ n< ^ a 
tight constraint on the model parameters a as < a < 1, 
where the other two model parameters are considered to 
500 K < T vir < 10 4 K and K w < 1/8. This result indicates 
that the high-redshift CSFH could be very flat. On the other 
hand, specifically, lower virial temperature T v i r leads to lower 
formation rates of Pop III stars, and thus the total formation 
rate is required to be higher in order to have more Pop I/II 
stars to produce emission. 

4 Summary and discussion 

In view of the inaccessible direct observations of high- 
redshift stars, especially Pop III stars and their hosts, we pro- 
pose an indirect method to study these high-redshift objects 
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by using the observed CIB emission, particularly, the NIRBE 
emission. In our model, the CSFH are described phenomeno- 
logically, in order to decrease the complexity of the model. 
Consequently, a strong constraint on the CSFH is found as 
< a < 1, if the virial temperature of the halos satisfies 
500 K < T Vl[ < 10 4 K. Such a flat CSFH basically agrees with 
the star formation rates inferred from high-redshift gamma - 
ray bursts up to redshift z ~ 8 - 9 |[T8l . In contrast, for higher 
redshift (e.g., z ~ 10), a recent measurement by 071 claimed 
that the star formation rate could become much lower. These 
observations imply that the realistic CSFH may exhibit a be- 
havior more complicated than a single power law. Anyway, 
on the other hand, the CSFH with < a < 1 predicts suf- 
ficiently bright sources to reionize the universe to account 
for the optical depth of cosmic microwave background pho- 
tons J38), which is measured to be r = 0.088 ± 0.015 by 
WMAP. 

The relationship between the parameters a and T YU can 
in principle be derived from the hierarchal formation model, 
although many uncertain astrophysical factors could be in- 
volved. Combining such a theoretical relationship with the 
result presented in this paper, we may obtain a more stringent 
constraint on the value of 7\,j r as well as the CSFH. 
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ence Foundation of China (grant nos 11103004, 11073008, 
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MOE of China (Grant Nos. CCNU12A01010). 
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